Correlation functions by Cluster Variation Method for Ising 
model with NN, NNN and Plaquette interactions 

E.N.M. Cirillo 

Universite Paris Sud, Batiment 425, 91405 Orsay Cedex, France. Emilio.Cirillo@math.u-psud.fr 

G. Gonnella, M. Troccoli 

Dipartimento di Fisica dell'Universitd di Bari and Istituto Nazionale di Fisica Nucleare, Sezione 
di Bari, via Amendola 173, 70126 Bari, Italy. E_mail: gonnella@ba.infn.it, troccoli@bari.infn.it 

A. Maritan 

International School for Advanced Studies, via Beirut 4, 34014 Trieste, Italy. 

Abstract 

We consider the procedure for calculating the pair correlation function in 
the context of the Cluster Variation Methods. As specific cases, we study the 
pair correlation function in the paramagnetic phase of the Ising model with 
nearest neighbors, next to the nearest neighbors and plaquette interactions in 
two and three dimensions. In presence of competing interactions, the so called 
disorder line separates in the paramagnetic phase a region where the corre- 
lation function has the usual exponential behavior from a region where the 
correlation has an oscillating exponentially damped behavior. In two dimen- 
sions, using the plaquette as the maximal cluster of the CVM approximation, 
we calculate the phase diagram and the disorder line for a case where a com- 
parison is possible with results known in literature for the eight- vertex model. 
In three dimensions, in the CVM cube approximation, we calculate the phase 
diagram and the disorder line in some cases of particular interest. The rele- 
vance of our results for experimental systems like mixtures of oil, water and 
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surfactant is also discussed. 
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I. INTRODUCTION 



Statistical systems are exactly solvable only in the simplest cases and often in unphysical 
spatial dimensions. In three dimensions, in particular, solutions of common statistical models 
like the Ising model are not known. On the other hand, numerical simulations sometimes 
are not a very efficient tool for studying systems with complex phase diagrams where an 
analytical control of the phase behavior of the system would be very helpful. If one is not 
interested in the critical behavior of the system, a very useful generalization of mean field 
methods for studying the phase diagram of spin systems is given by the cluster variation 



diagram that can be appreciated when a comparison is possible with exact or Monte Carlo 



properties of spin models. 

The CVM is based on the application of the variational principle of the statistical mechan- 
ics with the simplification that the density matrix of the system is factorized over clusters 
inside which the system is treated exactly. The size of these cluster, sometimes called "max- 
imal clusters" , depends on the accuracy requested and on the necessity of including all the 
interaction of the system and of reproducing the structure of the ground-states. Correlation 
functions have been calculated in the framework of CVM approximation only in few simple 
cases. Some properties of the pair correlation function of the two-dimensional Ising model 
have been studied in the pair (Bethe) and plaquette approximation levels of the CVmI. The 
CVM has been also applied to the study of the correlation function in the disordered phase 
of the two-dimensional ANNNI (Axial Next Nearest Neighbors Ising) mo dei. 

In this paper we present a procedure for the calculation of correlation functions that can 
be generally applied to any choice of maximal cluster in the CVM scheme. It can be also 
implemented in numerical calculations together with the Natural Iteration MethoJ which 
is a numerical minimization procedure very convenient for calculating the maximal cluster 
density matrix in the most complex cases. In particular we apply our procedure to study 



method 




This method generally gives a very accurate description of the phase 




It has been extensively used for the calculation of equilibrium thermodynamics 
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the behavior of the pair correlation functions in the paramagnetic phase of the Ising model 
with nearest neighbors (NN), next-to-the- nearest neighbors (NNN), and plaquette (P) inter- 
actions in two and three dimensions. In three dimensions, the phase diagram of this model 
can be correctly studied by implementing the CVM cube approximation with 8 independent 
site magnetizations in order to take into account the complex ground-state structure. Our 
goal will be to calculate the pair correlation function in the same approximation used for 
the phase diagram. 

The model with NN, NNN and P interactions has the same hamiltonian of the symmetric 
two-dimensional eight-vertex modeli. In three-dimensions this model is interesting also 
because it constitutes a realization of interacting random surfaces on a latticffl. Indeed, 
the Peierls interfaces between domains of different sign can be interpreted as an ensemble 
of interacting surfaces with a Boltzmann weight depending on the area, the mean curvature 
and the length of the intersection lines between the surfaces^. It has been shown0'0 that 
the phase diagram of this model, studied in mean field approximation, can well describe the 
phase diagram of experimental systems like fluid mixtures of water, oil and surfactant or 
other complex fluids. Recently, the model with NN, NNN and P interactions has been also 
put in relation with string theory and in particular with a discretized string model (the so- 
called gonihedric model) characterized by a zero energy cost for the area of surfaces^!. The 
CVM cube approximation has been applied to the study of the phase diagram of the model 
with NN, NNN and P interactions in the parameter region corresponding to the gonihedric 
mo del0. 

A question of particular interest in the phase diagram of models with competing inter- 
actions is the behavior of the pair correlation function in the paramagnetic phase. Here, 
generally, the so called "disorder line" separates a region where the pair correlation has 
the usual exponentially damped behavior from a region where the correlation develops an 
oscillatory exponentially damped behavioJil. In the case of the three-dimensional model 
with NN, NNN, P interactions, this latter region can be identified with the microemulsion 
phase of surfactant systems which is a disordered phase but with an ordered structure on 



short length-scaleal^. Therefore, in this model, the possibility of calculating the correlation 
functions in the same CVM approximation used for the phase diagram is interesting also for 
some physical applications. 

The plan of the paper is the following. In Sect. II we briefly summarize the CVM 
procedure and discuss how correlation functions can be calculated. In Sect. Ill we consider 
the bidimensional eight- vertex model. We calculate the disorder line and compare our results 
with other results existing in literature. In Sect. IV we will carry out the calculation of the 
correlation functions in the three-dimensional case and study the disorder line in two cases 
of particular physical interest. Some conclusions will follow. 

II. CVM AND CORRELATION FUNCTIONS 

In this section we first briefly describe the CVM approximation. We refer to the papers 
of refs.l for further details. Then we will discuss how the pair correlation functions can be 
evaluated starting from the CVM free energy. 

Suppose to have a spin hamiltonian defined on a given lattice as 

-/37^(a) = ^J«a« (1) 

where (3 is the inverse of the temperature, is the interaction parameter for the cluster a 
and cTq, is the product of spins cTj on the sites z G a. The first step of the CVM procedure is 
to choose a set of maximal clusters Ai in which the system is treated exactly such that each 
cluster a G X is included in one maximal cluster M E M.. Then the CVM approximated 
free energy of the system can be expressed! as the minimum of 

(3T = Jaia + a„Trp„ lnp„ (2) 

Here the first sum is the internal energy of the system with the being variational param- 
eters representing the expectation value of the multisite product of spins a a- The second 
term in the r.h.s of is the entropy of the system and is a sum over all the clusters which 
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are possible subclusters of the maximal clusters. The set of coefficients depends on the 
lattice structure and on the choice of the maximal cluster; they can be easily found applying 
the eqs.(16') in the paper by Ani. The trace Tr is the sum over the allowed configurations 
and Pa is the density matrix for the cluster a. In the case of Ising spins cTj = ±1 one can 
write 

Pa = 2-^-[l + Y.apif,] (3) 

/3Ca 

where Ha is the number of the sites in the cluster a, the sum extends over all the subclusters 
/3 of the cluster a, and ap = Hie/J'^j- ^ consequence of eq. (|) is that = liappp for a 
cluster (5 (Z a. The parameters ^q, have to verify the minimization conditions 

-f 

In order to study the pair correlation function for a group of spins o"^, it is convenient 
to introduce a local external field in the starting hamiltonian coupled to all clusters of the 
same type of 7. For one of these clusters the equilibrium condition becomes 

- f 

Therefore the connected pair correlation function is given by 

where we have used the linear response theorem and the equilibrium conditions (|^). Now, the 
free energy JF( J^, h^) can be considered formally as a function of the independent variables 
^-y instead of the h^. Therefore the last term of the above equation can be evaluated through 
its inverse matrix that can be calculated by differentiating (|^) with respect to 

dh d'^T 

Once that the above matrix is obtained, its Fourier transform can be simply inverted 
and will give the requested correlation function in Fourier space. However, except that in 
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the simplest cases, writing the matrix .4.^o,7r '^^^ ^ comphcate task from a computational 
point of view and requires some further considerations that we present here in a general 
way and that will become more transparent when applied to specific examples in the next 
sections. 

In the matrix ^70,7- there appear generally derivatives like d^a/d^-y (see, e.g., eqs. 
(|l^,|l^)) where 7 is a cluster of different type of a. The evaluation of these derivatives 



can be performed by differentiating the state equations for with respect to This 
procedure will give a system of linear equations 

in the variables d^a/d^^. Some of the equations (§) are homogeneous. More specifically these 
situations may occur: i) the cluster a belongs only to one maximal cluster. In such a case 
the equation (^) is homogeneous or not, depending if 7 is external or not to that maximal 
cluster, ii) The cluster a belongs to more than one maximal cluster and 7 is external to 
the maximal clusters with a in common. In such a case the equation (|^) is homogeneous, 
iii) The cluster a belongs to more than one maximal cluster and 7 belongs to one of the 
maximal clusters with a in common. In this case the equation is generally not homogeneous. 
It can be shown that a solution of the system @ can be obtained by setting equal to zero 
all the terms dC,a/dC,j appearing in homogeneous equations. This procedure will be more 
clear and it will be shown in explicit cases in the next two sections. 

Finally we observe that a great simplification occurs when we want to calculate the 
correlation function in the paramagnetic phase of a system with only even interactions. In 
this case only the derivatives d^a/d^-^, where a and 7 have both an odd (or even) number 
of sites, are different from zero. 

III. PAIR CORRELATION IN THE 8- VERTEX MODEL 



In the following of the paper our principal aim will be to calculate the two-spin correlation 
function and the disorder line of the model defined by the hamiltonian 
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<x,y> «x,y» 

where /3 is the inverse of the temperature and the ax are Ising variables defined on the sites 
of a cubic lattice. The three sums respectively refer to nearest, next-to-the nearest neighbors 
and plaquettes of the lattice. In this section we consider the CVM plaquette approximation 
of this hamiltonian on a square lattice and we use it to calculate the phase diagram, the 
correlation function in the disordered phase and the disorder line in the case J3 = 0. In the 
following section we will consider the three-dimensional case. 

3. 1 The plaquette approximation and the phase diagram. 

The natural choice for studying the phase diagram of the model in two dimensions is 
the CVM approximation level where the maximal cluster with independent density matrix 
is the plaquette cell of the square lattice. When a magnetic- field term is added to (||), the 
CVM free-energy density functional to be minimized is given by 

= ~ JlH<xy> Tr ((Ta^CTy p<a;y>) — ^2 E<<xj/>> "^"^{(^xf^y P«xy») + 

- JzY.-ai Tr((7xaj;(j^cr^pgnu,) - h,j:i{(Txpx) + E-n» Tr£(pjQu,)+ ^^^^ 

) + Ex Tr/:(p^.) 

where the spins in the argument of the traces are the spins in the cluster one is considering, 
C{a) = a log a for a real number a and pa is the density matrix related to the cluster of a 
given type a of the lattice; in the following we need also the matrix Pixyz] where [xyz] denotes 
a three site corner cluster. The density matrix p^ is subjected to the constraint Trp^ = 1 
and the smaller cluster density matrices are obtained by partial traces of p^. Notice that 
we have not assumed any a priori symmetry property for our density matrices; this implies 
that all possible states with different magnetic order in a single plaquette can be studied by 
this approximation. 

The plaquette density matrix p^ can be calculated by minimizing the free energy (|TUp 
via the natural iteration equationJ^; then, after partial traces, one can work out the phase 
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diagram of the system. In Fig. |T] the portion of the phase diagram with J3 = is shown. The 
hne separating the ferromagnetic and the paramagnetic phases is critical, while on the line 
separating the paramagnetic and the lamellar phase - sometimes called Super AntiFerromag- 
netic (SAF) phase - there is a tricritical point at Ji = 0.405 ±0.005, J2 = -0.4045 ± 0.0015. 
The existence of a first order phase transition between the paramagnetic and the lamellar 
phase is also confirmed by other authorsElS that have used the plaquette CVM approxi- 
mation to study the phase diagram of the 8-vertex model. Our results for the location of 
the tricritical point are in agreement with results of refs.0'Ei. The discrepancy of this result 
with the results of Monte Carlo simulations and perturbations methods where a nonuniver- 
sal critical behavior is found - for a list of references seelll, may be an artifact of the CVM 
approximation. 



3.2 The Correlation Function. 

For the calculation of the correlation function it is convenient to write the free energy 
(p!0|) as a function of the parameters 

= Tr (o-^.p^) 

Ixy — Tr((T3;Crj^p<a;j/>) 

Cxy = Tr(cr^(TyP«a;j^») (11) 
^xyz Tr ((Tj^CyCT^iPlxj/z] ) 

dxywz = T^T^i'^x^'y'^w'^zP^a^) 

that are related to the density matrices by: 
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1(1 + m^an, + m„cr^ + /^c^xcr^ 



(12) 



Px = 

P<xy> 4^-^ ' ""x'^x I '""y'-'y ' ''xy^x^yj 

P«xy» = ^(1 + f^x'^x + '^y'^y + CxyCxCy) 

P^n™ = ]^ (1 + "^^'^x' + "^i/CTy + m,^cr,^ + m^a^+ 

kyxw^y^x^w ~l~ ^xwz^ x^w^ z ~l~ ^wzy'^w'^ z'^y ~l~ ^zyx^z^y^x~^ 
dxywz(^x(^w(^ z^y^ 

We have seen in Sect. II that in order to calculate the pair correlation function, we need 
to evaluate the matrix 

Therefore we consider the state equation for the magnetization dJ^/drrii = 0: 



(13) 



f^^i = -jYl Tr(ailogp<,y>) + -^^Tr(ailogpi^n») + ^Tr(ailogpi 

<iy> iw 



(14) 



where the first and the second sums are respectively taken over the pairs of nearest neighbors 
and the plaquettes containing the site i. We start from the diagonal elements of the matrix 



(p!3|) given by 



dhj 

dlTLi 
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ETr 

<iy> 



,P<iy> 



+ 



+ — VTr 

256 



Pir-[ui 



my 



drrij 



where all the derivatives are evaluated in m. 



drrii 
\/z and 



(15) 



9mj 

has been assumed. 



Derivatives like ^ ^ 



„ , p , „ do not appear because they are odd functions of rrii evaluated 

ami ' ami ' ami ' 

d \/x, y, W, Z. 



at mj = 0. Moreover in the paramagnetic phase \xy — ^xy — ^ and dxywz 



To simplify the notation we introduce 

dkj 



"lyz 



drrij ' 



k,. 



rlh ■ 

^i^wiy 

drrii 
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and k'i 



dk 



yzw 



drrii 



(16) 



While the parameters /, c and d can be calculated as usual by solving the equilibrium 
state equations by means of the Natural Iteration Method, the calculation of ki, ^2 and 
is a little bit more intriguing. As explained in Sect. II it is convenient to consider the state 
equation obtained by setting equal to zero the derivative of the free energy with respect to 
h 



(17) 



Now, by differentiating this equation with respect to rrij with j = x,y,w, we obtain the 
system: 

dkx'y'z' 



Tr 



I yi—'z 



E 



ax'CTy/az' 



drrii 



j =x,y,w 



By solving this system of three equations in the variables ki, k2 and k^ one obtains 
A{C-D) , {C - D){B^ + BD -2A^) 



{D + B + 2A){D + B -2A) 



ko 



{D - B){D + B + 2A){D + B - 2A) 



(18) 



where 



2A\C + D-2B) + {B^ - DC){B + D) 
{D - B){D + B + 2A){D + B - 2A) 



(19) 



A = Tr ^ 



B = Tr 



C 



Tr 



D 



Tr 



Oxr-fW 



(20) 



The insertion of eqs. (|l^,^) in (|T5D gives the diagonal elements of (|T3D- Going further, 
by differentiating (|l^ with respect to irij with z,j nearest neighbours, one realizes that 



derivatives like 



with j ^ ^al appear. However, these derivatives are zero. Indeed, 



by differentiating the state equation for k^yz with respect to rrij with j external to the 
plaquette of [xyz], we obtain an equation similar to (|18D but without the first term aj in 



circular brackets. Changing the site j we again obtain a systems for the derivatives 



dhxy: 

dm-i 
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with j external to the plaquette of [xyz] , but in this case these derivatives are zero since the 
system is homogeneous. 

A consequence of the above considerations is that the matrix (ITBD can be written as 



with 



(3 



dhi 



(9m, 



7 = 1 



I -I' 



+ il + k2 + h 



7 ^ = y 

-7i <iy> 

-72 << iy » 
otherwise 

l + 2c + d 



;i + 4/ + 2c + d){l -Al + 2c + d) 



(21) 



1 2 

+ {l-k2 + h)- + (1 - fca' 



8kil 



l-2c + d 



1-d {l + 4:l + 2c + d){l-Al + 2c + d) 



(22) 



7i 



I 



21 



1-P {l + 4l + 2c + d){l-Al + 2c + d) 



+ ki 



l + 2c + d 



2(^2 + k3)l 



[l + Al + 2c + d){l-Al + 2c + d) 1-d^ 



1 + 4/ + 2c + d){l -4:l + 2c + d) 

(23) 



and 



_ 1 , . l + 2c + d 

- 4^^ + + '''>{l + 4l + 2c + d){l-4l + 2c + d) 



2kil 



4^ ^ ^^l + 2c + d 2(1 -d) (1 +4/ + 2c + d)(l -4/ + 2c + d) 



(24) 

The inverse of the Fourier transform of the matrix (|21|) gives the pair correlation function 
in Fourier space, or structure factor. It reads as 

1 



^(P) 



7-71 E;.=i cos(p/,) - 272(cos(pi + P2) + cos(pi - P2)) 



(25) 
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3.3 The disorder line. 

The hne where the coefficient of in the denominator of (p5| ) is zero is the socalled 
Lifshitz hne where the structure factor develops a maximum at a value of p different from 
zero. It is given by 71 + 272 = and is reported in Fig. |I[ The other interesting line for the 
behavior of the correlation function is the disorder line where the correlation in real space 
changes its behavior from a purely exponential decay to an oscillating exponentially decay. 
In order to calculate this line it is convenient to take the spherical average of the expansion 
of ( p5|) at small p so that the structure factor can be written as0: 

^^P) = 4 (26) 

1 + bp'^ + cp^ 

The spherical average can be simply realized by expanding (^) until the fourth power of p 
and taking Px = Py = p/ -\/20. The result is 

71 + 272 _ 1 71 + 872 ^27) 



7 - 4(71 + 72) ' 24 7 - 4(71 + 72) 

On the disorder line the zero of the denominator of (|25| ) change from pure imaginary to 
complex. This happens when 6^ — 4c = 0. In Fig. |I| we plot the disorder line obtained by 
eqs. (^. 

In some bidimensional model with competing interactions the disorder line coincides 
with a locus (One-Dimensional-Line) where the model can be solved exactly and has typical 
one-dimensional correlations^. The ODL line has been calculated for the model (§) with 
J3 = and is given byil 

p4J2 _|_ P-4J2 _|_ 9p-2J2 

2(1 + e^^.) 

It would be interesting to know whether this line, also reported in Fig. 1, coincides with 
the disorder line of the model (|^). In the paper^ the ODL line has been compared to the 
disorder line obtained in mean field approximation. Here we have an approximation which 
correctly reproduce the topology of the phase diagram of the model also at low temperatures. 
This allows us to compare our approximation for the disorder line with the expression (^Sf ) 
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at any value of Ji, J2- We see that the disorder hne calculated by CVM is very close to the 
ODL eq. (pHD, so that it is probably true that the two lines coincide. For further comparison 
we have plotted in Fig.l also the disorder line obtained by mean-field approximation. 



IV. THE 3D ISING MODEL WITH NN, NNN AND PLAQUETTE 

INTERACTIONS 

In this section we study in the CVM approximation the phase diagram and the disorder 
line of the model (P) defined on the cubic lattice. The ground states of the model @ 
are shown in Fig. 2. They can be simply obtained by minimizing for each value of the 
parameters Ji, J2, J?, the energy of a single cube0. When two or more cube configurations 
have the same minimum energy, it may be possible to construct a set of degenerate ground 
states by tiling the whole lattice with the degenerate cube configurations. We have chosen 
to study the particular cases J2 = J3 and J3 = — 2 J2 since in these cases the structure of the 
ground states is of particular interest for the applications0'0. The diagrams of the ground 
states for these cases are shown in Fig. 3. 

Due to the fact that the ground states can be expressed in terms of single cube configu- 
rations, it is necessary for a correct implementation of the CVM approximation to consider 
the cube with 8 independent site magnetizations as the maximal cluster in the CVM imple- 
mentation. Then the approximated free-energy is given by 

= ~ JlJ2<xy> Tr ((TxCTj^ P<xy>) ^ J2j2«xy» T^T^{o'xO'y p<:<:xy») + 

- -^^sEgn- Tr(cr^.o■J^c^^„cr^pJ□J.) - PJ2x hxTr^axPx) + Ec Tr£(pc) + 

(29) 

) - Ex Tr£(p^) 



where pc is the density matrix of an elementary cube. Also here, as in the bidimensional 
case, the coefficients of the entropic terms have been found by applying the prescriptions of 
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4-1 The phase diagram with J2 = J3 

This case has been studied in0'0 by using mean-field approximation and Monte Carlo 
simulations. In terms of dual surfaces, it corresponds to the case where the curvature, 
defined as the number of adjacent plaquettes forming a right angle, is not weig hted0. The 
phase diagram found by applying the natural iteration scheme to the CVM free-energy 
is shown in Fig. ^ At positive J2 the ferromagnetic phase is stable. On the transition line 
separating the ferromagnetic and the paramagnetic phases there is a tricritical point at J**" = 
0.029, J2 = 0.0635. The coordinates of the tricritical point in mean field approximation 
and by Monte Carlo simulations are respectively given by J^{'^^^ = 0.11, j!^'^^^ = 0.0275, 
jtr,MC _ g J^^'^^^' = 0.064. A relevant difference between the mean-field and the CVM 
phase diagram is that in the CVM approximation the paramagnetic phase extends until 
zero temperature, in agreement with results of Monte Carlo simulations. At negative J2 
the ordered stable configurations can be constructed starting by the cubes in Fig. |^. The 
phase between the phases 4 and 7 in Fig. 4 has the same configuration of phase 7 but 
with the magnetizations m = 0. This phase is separated from the phases 4 and 7 by a 
first-order and a critical line, respectively. The phases 4 and 7 are bicontinuous in the sense 
that the domains of spins of one sign form a connected network invading all the lattice and 
intertwined with the network formed by the spins of the other sign0. Also mean field and 
Monte Carlo simulations show stable bicontinuous configuration in this region. However, 
there are discrepancies between the phase diagrams obtained by different methods at J2 < 0. 
The origin of the discrepancies is probably in the particular nature of the ground states at 
J2 < — I Jil/4. As it can be seen from Fig. |^, the cube configurations 4,4, 7 are degenerate 
for J2 < — I Jil/4 and an infinite number of bicontinuous ground states can be built up using 
these configurations. Here, probably, a low-temperature expansion is needed to understand 
the correct nature of the stable phases, but this will be the matter for a future study. It is 
interesting to observe that on the transition line between the paramagnetic phase and the 
4 phase at negative J2 there is a tricritical point located at Ji = 0.87, J2 = —0.32. It would 
be interesting to check the existence of this tricritical point by Monte Carlo simulations. 
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4-2 The disorder line. 

In three- dimensions the calculation of the disorder line proceeds in the same two steps 
as in the bidimensional case but with complications due to the higher number of variational 
parameters in the pc expansion. This implies the existence of many terms like the quantities 
ki , k2 , etc. of the previous section to evaluate. 



The first step is the calculation of the elements of the inverse correlation matrix (|13|), 
that in tree dimensions becomes 



P 



dhi 



dm. 



y / m,=oyz 



1 i = y 
-71 <iy> 
\ -72 << ly » 

(30) 

-73 «< iy »> 
otherwise 

The explicit expressions of the coefficients 7^ are given in Appendix A for clearness purposes. 
Going on with the second part of the calculation, we use these matrix elements to compute 



the coefficients (|27|) of the small-p expansion of the structure factor, that in three dimensions 
are given by: 

71+ 472 + 473 1 7i + I672 + 2873 



(31) 



7 — 671 — 1272 — 873 36 7 — 671 — 1272 — 873 

Also here, as in the bidimensional spherical average has been performed taking 

Px = Py = Pz = pI "\/3- When 6^ — 4c < the correlation function in real space is given by 



G(r) 



const _„/f . 27rr 
e '''sm— ^ 



(32) 



where 



5 



27r 



^ 4c 



1/2 



_i b_ 

2^c 4c 



1/2 • 



(33) 



The disorder line corresponds to the condition 6^ — 4c = 0. The results for the disorder line 
and for the Lifshitz line 6 = are shown in Fig. ^. At low temperatures both the lines 
behave as J2 = — Ji/4. 
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4-3 The phase diagram at J3 = — 2J2. 

The ground-state structure corresponding to this case is shown in Fig. 3. We see that the 
ferromagnetic phase is bounded from a phase with spin of different sign in alternate planes. 
We call this last phase a lamellar phase. At finite temperature the phase diagram has been 
studied in mean field approximation in a previous paperS. It can well describe the phase 
diagram of experimental mixtures of oil, water and surfactant. The CVM approximation 
is the same used before and the calculation of the disorder line also proceeds in the way 
already discussed in the preceding paragraph of this section. Therefore we can go on by 
describing the phase diagram shown in Fig. |^. 

At low temperatures and positive J2 the stable phase is the ferromagnetic one, while 
for negative J2 we find the lamellar phase being stable. On the transition line between 
the ferromagnetic and the paramagnetic phases there is a tricritical point at Ji = 0.273, 
J2 = —0.049, with the first-order part of the line joining the boundary of the lamellar 
phase region at Ji = 0.305, J2 = —0.077 (4-phase point). This feature, together with the 
observation that the disorder line intersects the first-order ferro-paramagnetic transition 
line in Ji = 0.29, J2 = —0.065, amidst the tricritical and the 4-phase point, well represents 
the experimental fact that the microemulsion phase (corresponding to the region of the 
paramagnetic phase below the disorder line) can coexist with the ordered homogeneous 
phase^3@. We also observe that the Lifshitz line intersects the ordered phases at the point 
Ji = 0.2965 ±.0005; J2 = -0.00705 ± .00005. The region close to the 4-phase point is shown 
in the inset of Fig. 7. 

V. CONCLUSIONS. 

The CVM approximation is generally used for studying the topology of phase diagrams 
of spin systems. In this paper we have focused on the calculation of the correlation functions 
in the framework of this approximation. This calculation can become complicate expecially 
in the case of three-dimensional systems and we have discussed the origin of these com- 
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plications. As particular models we have considered the Ising model with NN, NNN, and 
plaquette interactions on the square and on the cubic lattice. We have calculated the pair 
correlation function in the paramagnetic phase of these models. We have also given new 
results about the phase diagrams of these models for choices of parameters useful for the 
description of experimental systems. 

In particular, in the two-dimensional case with zero plaquette interaction, we have cal- 
culated the disorder line in the paramagnetic phase and we have seen that it is very close to 
the One-Dimensional Line where the model is exactly solvable. This confirms the conjecture 
that the disorder line and the ODL line coincide also in this model. In the three-dimensional 
case we have chosen two planes in the parameter space where the phase diagram may have 
a relevance for the description of systems of interacting surfaces with fluctuating topology. 
Fluid mixtures of oil, water and surfactant are an example of these systems. We have cal- 
culated the disorder line and the Lifshitz line which limit regions of the paramagnetic phase 
detectable through scattering experiments. Our results, shown in Fig. confirm that the 
three-dimensional Ising model with NN, NNN, and plaquette interactions can describe many 
of the experimental features appearing in complex fluid mixtures. 



VI. APPENDIX A 



The matrix (|30D can be evaluated through the same procedure shown in sect. Ill, starting 
from the single spin magnetization state equation ([T^), that in 3D turns out to be: 



(^hi = (f^i log Pc)-^Y1 Tr(f^i log p^nj 



c9« s n™ 



+ 7 Tr(ailogp<ij^>) - ^Tr(ailogpi) (34) 

^ <iy> ^ 

Differentiating this equation with respect to the site magnetizations, we realize that we have 
to deal with a long list of derivatives of the kind: = where ^q, = Tr((TaPa) is the 
parameter in the expansion of the density matrix corresponding to the subcluster a. By 
symmetry considerations and taking into account that in the paramagnetic phase the only 



derivatives different from zero are those from subclusters with an odd number of sites, we 
can reduce the number of derivatives to 33, hsted in Fig. ||. Repeating the procedure used 
in section III to obtain (0), we can get to a system of 33 hnear equations which, although 
very long, can be straightforwardly solved and it will not be reported here. The solution of 
this system gives us the coefficients fcj , , bi appearing in Fig. g and allows us to calculate 
the expression of the matrix elements (^), that we write down here explicitly: 

^ = 8 [ + 6hC^ + 3k2C^ + 3hC^ + 3hC^ + 3kr,C^ + 3keC^ + 3/^7^^ + 

+6ksC^ + 6kgC~^ + kioC^ + 3knC^ + 3kuC~^ + k^C^ + GaiC^'' + 3a2C^^ + 
+3a3C^ + eauC^ + 3ai5C^ + 3aieC~^ + 304^^" + 3a5C^' + 3aeC^'' + 3a7C~^ + 
+6a8C^ + 6agC~^ + aioC=^' + 3auC^'' + 3ai2C~^ + aigC" + + 362^^" + 
+363^^' + &4C^' ] - 12 [ + 2kiP^ + k2P^ + A;3P° ] + 6^°^= - 1 (35) 

- 7i = 4 [C^ + 2ki{C^ + C^ + C^) + k2{2C^ + C~^) + k^(2C'^ + C^) + k^^C^ + 2C^) + 

+k^{C^ + 2C^) + kQ{2C^ + C~^) + k^^C"^ + 2^3) + 2ks{C^ + + C~^) + 2kg{C^ + 

+ C^) + fcioC^ + A;n(C^ + C^) + fci2(2C2 + C°) + fcisC^ + 2ai(C2^ + + C^) + 
+a2(2C3^ + C^) + a3(2C3^ + C^) + 2014(6'^'= + + C^) + ai^{C^' + 2C') + 
+ai6(C6 + 2C^) + a4(C^ + 2C^) + 05(6-' + 2C') + a6(2C^ + C^) + a^{C'' + 2C^) + 
+2a^{C^ + C^ + C^) + 209(6-=^ + C' + C^) + aioC^ + aii(C^ + C^) + 
+ai2(2C2 + C°) + ai3C^ + hiC^' + 62(C° + 26'^'=) + 63(2C=^' + C^") + WC^''] + 
- 4 [ + A;i (P2 + P°) + A;2P^ + A:3P^ ] + L° (36) 

- 72 = 2 [C^ + 2fci(C=^ + + C^) + fc2(C6 + 20') + A;3(C2 + 2C^) + A;4(2C^ + C^) + 
+k^{C^ + 2C5) + A;6(C^ + 2C^) + A;7(C2 + 2C') + 2A;8(C2 + + C^) + 
+2kg{C^ + C^ + C^) + A;ioC^ + A;n(2C2 + C°) + fci2(C^ + 2C^) + k^C^ + 
2ai(C3^ + + C^) + a2(C6 + 2^7) + a^{C^^ + 2C^) + 2ai4(C2" + + C^) + 
+ai5(2C=^^ + C^) + ai6(2C3'= + C^) + a^{2C' + C^) + 05(^2^ + 2C^) + 
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aioC'^ + aii{2C^' + C°) + a^2{C^ + 2C^) + aiaC'^ + feiC'^ + 62 (C^^ + 2C=^^) + 
+63(C° + 2C=^^) + 64C^^] - + 2A;iPi + ^2^° + fcsP^] (37) 

- ^3 = + 6A;iC^ + SfcaC^ + ?,k^C^ + 3/^4^^°' + SfcgC^ + S/tgC^^ + 3/^7^^ + Qk^C^ + 
+6fc9C^ + A;ioC° + 3A;iiC^ + +3/^12^^ + k^^C^ + 6aiC^ + 802^3 + 303^^ + 
+6ai4C^" + SaisC^ + SaigC^" + 804^°" + SagC^ + SagC^^ + ^a^C^^ + GagC^ + 
+6a9C=^ + aioC° + SanC^ + 3012^^" + aigC^ + + S^aC^^ + SfogC^" + 64^° 

(38) 

where we used the notations: P"^ = ^ pi _ 1 p2 _ 

256 \^ y 256 Y'^p^ ) 

256 y'^^D'^ y \/'5n^/' \P<a:y>/' 16 VP<i!H>/ 

C° = (2^y2 Tr (^^^ , a being the name of the cluster as they are classified in Fig. 
and a a = Iljea in Section II. 
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FIG. 7. 
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Figure captions 



Fig.l: Phase diagram of model (§) in D = 2 with J3 = 0. The dashed lines are critical lines 
separating the disordered phase from the ferromagnetic and the lamellar (SAP) phases. The 
first-order transition between the paramagnetic and SAP phases is represented by a solid 
line. The dotted line is the one-dimensional line of equation (pH]), while the light dashed and 



dash-dotted lines are the mean field and plaquette approximation calculation of the disorder 
line, respectively. 

Fig. 2: Ground states of the 3D model. The dot • indicates a = +1 while a = —1 in the 
remaining sites. The labels refer both to the phase and to the dotted subclusters. The 
complementary (non-dotted) subclusters can be referred to adding a "c" after the label. 

Fig. 3: Ground states of the 3D model for the values J2 = J3 and J2 = —^3/2 of the coupling 
constants. 

Fig. 4: Phase diagram of model in D = 3 with J3 = J2. The dash-dotted line is the 
disorder line in cube approximation. The dotted line is the Lifshitz line. In the narrow strip 
between the 7 and 4 phases we find a type-7 phase with the magnetizations m of Pig. 5 
equal to zero. As before, solid and dashed lines represent first-order and critical transition 
lines. 

Fig. 5: Phases at finite temperature appearing in the diagram of Pig. ^. Capital letters 
indicate bigger magnetizations. 

Fig. 6: Derivatives = |^ with respect to the magnetization mj in the site circled (see 
App. A for the notations), ki, ai and bi are respectively the derivatives of the 3-, 5- and 
7-sites subclusters, represented by the dots. /^Prom ki{ai) to ^3(03) the cluster considered 
is a = 5(5c); from ^4(04) to kgi^ag), a = 5(5c), from fcio(aio) to A;i3(ai3), a = i(ic), from au 
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to ai6, a = 5c and there is no correspondent in the fcj's; from bi to 64, a = Ic. 

Fig. 7: Phase diagram of model in D = 3 with J3 = — 2J2. The dash-dotted hne is the 
disorder hne; the dotted hne is the Lifshitz hne. In the inset, the region where the disorder 
and the Lifshitz hnes cross the first-order Ferro-Paramagnetic transition hne is enlarged. 
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